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摘 要 : 本 研究 基于 CASA 模 型 获得 黄土 高 原 2000 一 2018 年 植被 NPP 数 据 集 ,通过 相关 性 ` 残 差分 析 等 方法 ,定量 
分 析 了 和 气候 变化 和 人 类 活动 对 植被 NPP 变 化 的 影响 。 结 果 表 明 :(1) 黄土 高 原 2000 一 2018 年 植被 NPP 呈 东南 高 
西北 低 分 布 格局 ,研究 区 约 86.86% 区 域 植被 NPP 时 增加 趋势 ,主要 为 退耕 还 林 还 草 的 核心 区 ;时 下 降 趋 势 的 面积 占 
13.14% , 主要 分 布 在 西北 部 的 干旱 区 。(2) 气候 变化 和 人 类 活动 对 植被 NPP 变 化 的 贡献 率 分 别 为 48.78% 和 51.22%， 
在 空间 上 有 明显 的 差异 性 , 人口 稀 少 、 发 展 相对 落后 地 区 的 植被 变化 主要 受气 候 变 化 的 影响 ;人 口 众多 、 经 济 发 展 
较 好 地 区 植被 变化 以 人 类 活动 为 主 。(3) 黄土 高 原作 为 干旱 半 干 旱 区 植被 对 气候 变化 的 响应 较为 敏感 , 随 着 人 类 活 
动 的 日 益 频 繁 ,气候 和 人 类 活动 共同 作用 于 黄土 高 原 的 植被 变化 。 本 研究 有 助 于 理解 气候 变化 和 人 类 活动 对 植被 
动态 变化 的 影响 ,并 为 黄土 高 原 的 植被 恢复 和 高 质量 发 展 提供 科学 依据 。 
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植被 是 陆地 生态 系统 的 主要 部 分 ,是 连接 大 
气 \ 土 壤 、 水 分 的 自然 纽带 ,是 反映 生态 环境 的 重要 
“指示 带 ”"”。 植 被 净 初 级 生产 力 (Net Primary Pro- 
ductivity,NPP) 是 植被 在 单位 时 间 .单位 面积 上 通过 
光合 作用 累积 的 有 机 物 的 数量 ” ,是 陆地 生态 系统 
中 碳 循环 的 重要 部 分 ,可 以 直接 反映 全 球 或 区 域 生 


究 。 在 全 球 变 暖 的 大 背景 下 ,气候 变化 对 植被 NPP 
的 影响 一 直 是 研究 热点 ;同时 ,社会 政策 变化 、 
城市 扩张 ”等 人 类 活动 的 作用 也 是 影响 植被 变化 
的 重要 因素 之 一 ,因此 ,探究 气候 变化 和 人 类 活动 
对 植被 动态 的 影响 ,在 未 来 各 区 域 应 对 气候 变化 和 
国家 生态 恢复 政策 等 方面 具有 积极 意义 。 


态 系统 的 质量 状况 ?9 ,在 生态 系统 碳 源 、 碳 汇 的 调 
节 过 程 中 发 挥 重要 作用 ”*。 因 此 ,对 植被 NPP 时 空 
变化 和 影响 因素 的 探究 ,不 仅 可 以 对 植被 的 变化 和 
了 驱动 进行 深入 的 理解 ,也 能 为 当地 生态 环境 保护 和 
可 持续 发 展 提供 科学 依据 。 

20 世 纪 60 年 代 , 主要 通过 实测 获取 植被 NPP 数 
据 , 但 这 种 方法 效率 低 且 难以 在 大 尺度 上 推广 ”9 。 
之 后 随 着 遥感 技术 的 发 展 , 通 过 遥感 模型 估算 长 时 
间 序 列 植被 NPP 的 方法 推广 开 来 汪 ,国外 学 者 主要 
基于 气候 统计 模型 2 生态 机 理 过程 模 型 所 和 遥感 
光 能 利用 率 模型 “开展 对 植被 NPP 的 研究 。21 世 
纪 之 后 , 光 能 利用 率 模 型 (CASA) 被 广泛 应 用 于 各 种 
JR BE FE BK NPP FY zs 2E [bU 7 OK oh R BK) O 
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黄土 高 原生 态 环境 脆弱 ,是 黄河 流域 泥 沙 的 主 

要 来 源 地 ,也 是 我 国 水 土 保持 重点 保护 区 ,其 生态 
系统 的 变化 关系 ee er 
生态 屏障 的 工程 效益 。 为 改善 黄土 高 原 日 益 恶 化 
的 生态 环境 ,从 1999 年 开始 实 ae 
RIESE TRE ,植被 变化 及 其 驱动 因素 是 
该 阶段 研究 的 重要 内 容 。 一 方面 ,侧重 于 气候 因素 
对 植被 变化 的 影响 , 韩 丹 丹 等 ”研究 黄土 高 原 植被 
变化 对 极端 气候 的 响应 ,结果 表明 NDVI 变 化 与 极 
气温 指数 存在 相关 性 ,而 与 极端 降水 指数 未 表现 
出 明显 的 相关 性 ; 程 晓 奢 等 ”研究 生态 恢复 与 局 地 
气候 变化 关系 ,结果 显示 黄土 高 原 气 候 受 全 球 气候 
变化 的 影响 ,植被 覆盖 度 的 变化 也 会 影响 局 地 温 
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BE .风速 .湿度 的 变化 。 此 外 ,植被 变化 还 与 人 类 活 
动 存 在 密切 关系 ,城镇 化 发 展 扩 .生态 恢复 工程 的 
实施 宝 同 样 影响 着 植被 变化 。 但 是 当前 研究 多 基 
于 NDVI 等 植被 指数 ,而 植被 NPP 能 直接 反映 植被 
光合 作用 积累 的 有 机 物 数量 ,对 气候 和 人 类 活动 变 
化 更 为 敏感 , 且 目 前 对 人 类 活动 的 影响 研究 多 基于 
定性 .比较 的 方法 ,缺少 对 人 类 活动 贡献 率 的 定量 
研究 。 因 此 ,本 文 利用 CASA 模型 计算 2000—2018 
年 植被 NPP, 基 于 相关 性 分 析 、 残 差分 析 等 方法 , 定 
量 探究 气候 变化 和 人 类 活动 在 植被 NPP 变 化 中 的 
相对 贡献 率 , 以 期 为 黄土 高 原生 态 保 护 策略 的 制定 
和 黄河 流域 的 高 质量 发 展 提供 科学 依据 。 


1 研究 区 概况 与 方法 


1.1 研究 区 概况 

3& E BR (107928 ~111°15’E, 35°21’ ~39°34'N) 
位 于 我 国 中 部 地 区 ,包含 内 蒙古 .宁夏 .陕西 .山西 、 
甘肃 河南 和 青海 7 个 省 (自治 区 ) ,面积 约 为 64x10* 
km( 图 1)。 地 势 西 高 东 低 ,属于 干旱 半 干 旱地 区 ， 
年 均 温度 为 4~14 % ,降水 量 为 200~800 mm。 由 于 
干旱 少雨 的 气候 、 沟 蜜 纵横 的 地 貌 SY AE JA 
类 活动 的 干扰 等 原因 ,导致 黄土 高 原 成 为 水 土 流失 
最 严重 的 地 区 和 黄河 流域 的 主要 来 源 地 。 为 了 改 
善 黄土 高 原 的 生态 环境 ,在 1999 年 开始 大 规模 的 实 
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经 有 了 较 大 的 改善 。 
1.2 数据 来 源 及 处 理 

采用 朱文 泉 等 “改进 的 CASA 模型 ,利用 其 在 
IDL 平 台 下 开发 的 模型 插件 ,实现 对 研究 区 植被 
NPP 的 估算 。 模 型 佑 算 所 需 数据 包括 :(1) 归 一 化 
植被 指数 (NDVD) 来 自 中 国 科 学 院 资 源 环境 科学 与 
数据 中 心 (https://www.resdc.cn/) 的 月 尺度 植被 指数 
空间 分 布 (1998 一 2020 年 ) 数 据 集 。(2) 植被 类 型 数 
据 采用 GLC2000 项 目 开发 的 SPOT 4 全 球 土地 覆盖 
数据 中 国 子 集 (https://data.tpdc.ac.cn/zh-hans/)。(3) 
气象 数据 选取 2000 一 2018 年 月 平均 气温 和 月 总 降 
水 量 站 点 数据 来 自 中 国 气 象 数据 网 (http://data.cema. 
cn/) ;太阳 总 辐射 数据 采用 国家 青藏 高 原 科 学 数据 
中 心 (https://data.tpdc.ac.cn/zh-hans/) 的 中 国 区 域 地 
面 气 象 要 素 驱 动 数据 集 (1979 一 2018 年 ), 利 用 Mat- 
lab 工具 进行 格式 转换 FU TEER SH o 
(4) DEM 数据 从 地 理 空间 数据 云 获 取 (http: /www- 
gscloud. cn)。 对 上 述 所 有 数据 进行 投影 转换 ER 
样 裁剪 等 预 处 理 , 均 统一 为 1 kmx1 km 分 辨 率 、 
WGS_1984_Albers #252. 


2 研究 方法 


2.1 CASA 模型 
植被 NPP 作 为 生态 系统 服务 功能 一 个 重要 指 


施 退 耕 还 林 还 草 `. 治 沟 造 地 等 生态 工程 政策 , 随 着 
生态 保护 工程 的 实施 ,目前 黄土 高 原 的 生态 环境 已 
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图 1 研究 区 示意 图 


Fig. 1 Location ofthe Loess Plateau 


标 , 反 映 了 气候 环境 以 及 人 类 活动 对 植被 的 综合 
作用 ,CASA 模 型 由 于 参数 易 获 取 、 精 度 较 高 过程 
简单 等 优点 ,已 被 广泛 用 于 植被 NPP 的 计算 
NPP(«,t) = APAR(x, t) X e(x,t) (1) 
APAR(x,t) = SOL(x,t) X 0.5 X FPAR(x, t) (2) 
&(x,t) =T (x,t) X T (nt) x W (x,t) Xe, (3) 
式 中 ;NPP(w) 为 像 元 x 在 1 月 份 的 植被 净 初 级 生产 
力 ;APAR(x) 为 像 元 x 在 1 月 份 获得 的 光 和 有 效 辐 
射 ;SOLG.0 为 像 元 x 在 :月 份 的 太阳 总 辐射 量 ;FPAR 
(x,t) fot TE t H BOAT ARDEA CHAT AY 
吸收 比例 ;e(x) 为 植物 把 吸收 的 光合 有 效 辐 射 转化 
为 有 机 碳 的 效率 ;T(x) 为 在 低温 和 高 温 环 境 下 植 
物 自 号 生理 作用 对 光合 作用 产生 限制 而 降低 其 净 
初级 生产 力 的 程度 ; T(x) 为 植物 在 生长 最 适 温度 
向 高 温和 低温 转变 时 光 能 利用 率 逐 渐变 小 的 趋势 ; 
Wt) AK A RSS GRE RU SC B BR le RE ses. 
为 理想 状态 下 最 大 光 能 利用 率 。 
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2.2 年 际 变化 趋势 分 析 

采用 一 元 线性 回归 方程 计算 植被 NPP 的 年 际 
变化 趋势 ,并 将 其 斜率 定义 为 植被 NPP 的 年 际 变 
化 趋势 率 (slope) ,公式 如 下 : 
nx V (ixNPP)) - Yi NPP, 

nx 2n - 2i 

式 中 :slope 为 NPP 与 时 间 变 量 拟 合 的 一 元 线性 回归 
方程 的 斜率 ;i 为 时 间 变 量 , 等 于 1 到 nn 的 整数 ;n 为 
研究 时 段 的 年 数 ,等 于 19;NPP, 为 第 ;年 的 生长 季 平 
EJ NPP, slope<O 和 slope>0 分 别 表示 NPP 序列 随时 
间 递 减 和 递增 ;slope 绝对 值 越 大 ,表明 NPP 变化 越 
快 。 变 化 趋势 的 显著 性 采用 检验 ,公式 如 下 : 


slope = (4) 


F-Uxn-2 5 
x 0 (5) 
U= (5-9) (6) 
Q= 3.-5) 7) 


SUH EUN VE2EOE TT A; HPA All; y BEAN i 
年 NPP 的 实际 值 ; 了 表示 第 i 年 的 回归 值 ; y 表示 多 
年 NPP 的 平均 值 ;n 表 示 时 间 序 列 。 
2.3 相关 分 析 

相关 分 析 是 揭示 两 个 及 以 上 处 于 同等 地 位 的 
随机 变量 间 关 系 的 分 析 方 法 ,本 研究 基于 像 元 分 析 
分 别 计算 NPP 与 不 同 气候 因子 之 间 的 关系 ,相关 系 
数 的 计算 如 下 


n 


Ys. — ky; =¥) 


R= = (8) 


n 


Èc- Èo- 
SOAR, AM KAM n AEA S x yal xy AY 
第 i 个 值 ; x y 分 别 是 变量 x、y 的 平均 值 。 相 关系 
数 介 于 -1~1 之 间 , 其 绝对 值 越 接 近 1, 两 变量 的 相 
关 程 度 越 强 , 越 接近 0, 相 关 程 度 越 弱 。 
2.4 残 差 分 析 

在 数理 统计 中 , 残 差 分 析 是 指 实际 真 值 与 预测 
值 之 间 的 差 值 ;在 生态 学 中 , 残 差 分 析 是 通过 建立 
回归 模型 ,具体 量化 气候 变化 和 人 类 活动 对 植被 变 
化 的 影响 后 。 主 要 分 两 个 步 又:(1) 以 降水 .气温 时 
间 序 列 数据 为 自 变 量 ,基于 回归 模型 拟 合 得 到 的 
NPP 预测 值 , 即 NPPce( 代 表 气 候 变 化 对 植被 NPP 的 


影响 );(2) 实 际 真 值 NPP, 和 NPPe 的 差 值 为 NPP 残 
2$ , Bl NPP (代表 人 类 活动 对 植被 NPP 的 影响 ) 公 
式 如 下 : 

NPP, =aXT+bxPte (9) 

NPP,,, = NPP p - NPP oc (10) 
式 中 ;NPPcc 和 NPP, 分 别 表示 基于 回归 模型 的 NPP 
预测 值 和 基于 遥感 影像 的 NPP 观测 值 ;NPP 为 残 
差 ;T 和 PP 分 别 表 示 平 均 气温 和 降水 量 ;a、b\e 为 模 
型 参量 。 当 NPPns > 0 表示 人 类 活动 对 NPP 有 积极 
影响 ; 当 NPPm=0 表 示人 类 活动 对 NPP 基 本 没有 影 
Wf] ; 24 NPPu. <0 表 示人 类 活动 对 NPP 变化 有 消极 影 
响 。 由 于 本 文 不 对 其 他 因子 的 贡献 率 做 分 析 , 即 认 


为 气候 变化 贡献 率 + 人 类 活动 贡献 率 =1。 
3 结果 与 分 析 


3.1 黄土 高 原 植被 净 初 级 生产 力 时 空 变化 

2000—2018 年 黄土 高 原 植被 NPP 呈 现 南 高 北 
低 、 东 高 西 低 的 分 布 格局 (图 2a) ,其 中 ,高 值 区 主要 
分 布 在 南部 的 关中 地 区 和 东部 的 山西 省 ,2012 年 达 
到 最 高 值 ; 低 值 区 主要 分 布 在 西北 部 的 内 蒙古 自治 
区 、 陕 西 省 北部 .宁夏 回族 自治 区 和 甘肃 省 部 分 地 
区 。 由 图 2b 可知 ,植被 NPP 均 值 年 际 变化 较 大 ， 
2000 年 和 2001 年 均 小 于 350 g-m” - a^ ,其 余年 份 的 
值 均 处 于 400~500 gem? a”, EII EFRA , 增 
长 速率 为 6.32g'm”.a 。 将 黄土 高原 植被 NPP 均 值 
所 占 面积 百分比 划分 为 4 个 等 级 :0~250 gema, 
250-500 g* m^*a,500-750 g*m^-a .>750 gm a'。 
结果 表明 ,植被 NPP 均 值 在 0~250 e+ m+ a 的 面积 
占 比 在 2000 年 和 2001 年 分 别 为 40.74% 42.73% ,其 
余年 份 均 小 于 35% , 且 在 整个 研究 时 限 内 该 区 间 面 
积 占 比 有 减少 的 趋势 ;植被 NPP 均 值 在 250~750 g- 
ma 区 间 的 面积 占 比 最 大 约 为 53%~77% ;植被 
NPP 均 值 在 >750 g- m^ - a 的 面积 占 比 2000 年 和 
2001 年 小 于 5% ,其 余年 份 面积 占 比 在 6%~14% 波 动 
变化 , 且 呈 缓慢 上 升 的 态势 。 整 体 上 ,2000 一 2018 
年 黄土 高 原 植 被 NPP 呈现 西北 低 东 南 高 的 分 布 格 
局 和 波动 增长 的 变化 趋势 。 

2000—2018 年 黄土 高 原 植被 NPP 拟 合 斜率 在 
-39.44~36.55 g-m^-a (图 3), 呈 增加 趋势 的 面积 
比 为 86.86% ,分 布 于 研究 区 的 大 部 分 地 区 ,其 中 , 极 
显著 增加 和 较 显著 增加 的 区 域 主要 分 布 在 陕西 北 
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(b) NPP 均 值 等 级 /(g:m?:a 1!) 
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图 2 黄土 高 原 2000 一 2018 年 植被 NPP 均 值 (a) 及 各 区 间 面 积 百分比 (b) 
Fig.2 Changes in annual NPP (a) and percentage of area in each interval (b) in the Loess Plateau from 2000 to 2018 
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(b) 显著 性 检验 
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图 3 黄土 高 原 2000 一 2018 年 植被 NPP 变 化 趋势 (a) 及 显著 性 检验 (b) 
Fig. 3 Trends of NPP (a) and significance test (b) in Loess Plateau from 2000 to 2018 


部 (延安 市 .榆林 市 ) .山西 吕梁 山区 一 带 ( 朔 州 市 、 
忻州 市 .吕梁 市 太原 市 ) 宁 夏 的 固原 市 和 甘肃 东 
南部 (临夏 市 定西 市 天 水 市 . 庆 阳 市 ); 变 化 不 显 
著 的 区 域 主要 分 布 在 内 蒙古 (鄂尔多斯 市 . 呼 和 浩 
特 市 5 25 p 8.) .山西 省 晋 中 市 .宁夏 .甘肃 ( 白 
银 市 .兰州 市 ) 和 青海 省 ; 呈 减 少 趋势 的 面积 占 比 为 
13.14% , 主要 零星 分 布 于 黄土 高 原 东 南部 的 河南 
省 .陕西 省 关中 地 区 和 西北 部 的 边远 地 区 。 
3.2 气候 变化 对 植被 净 初 级 生产 力 的 影响 

IPCC 评 佑 报告 指出 ,全球 气候 正在 变 暖 ,而 且 
气候 变 暖 与 碳 循环 存在 显著 的 正 相 关 关 系 ,但 是 不 
同 地 区 气候 因子 对 植被 NPP 积 累 的 作用 不 尽 相 
同 。 由 图 4a 可知 ,黄土 高 原 的 累计 降水 量 与 年 均 植 
被 NPP 的 相关 系数 为 0.47, 其 中 正 相关 关系 占 整 个 


研究 区 的 93.79% , 主要 分 布 在 黄土 高 原 北 部 的 内 蒙 
古 自 治 区 BEJE .山西 .甘肃 .青海 和 宁夏 回族 自治 
区 的 大 部 分 地 区 ; 负 相 关 关 系 占 整个 研究 区 的 
6.21% , 主要 分 布 在 黄土 高 原 南部 的 关中 平原 .河南 
省 。 由 图 4 可知, 黄土 高 原 的 年 平均 气温 与 年 均 植 
被 NPP 的 相关 系数 为 -0.05 ,其 中 正 相 关 关 系 占 整个 
人 研究 区 的 41.5% , 主要 分 布 在 陕西 榆林 .山西 西部 
CB Xe 太原 市 . 晋 中 市 和 忻州 市 ) 宁夏 和 青海 
省 ; 负 相 关 关 系 占 整个 研究 区 的 58.5% , 主要 分 布 在 
内 蒙古 自治 区 鄂尔多斯 市 .陕西 省 (延安 市 .渭南 
市 ) 和 甘肃 省 定西 市 。 总 体 上 ,黄土 高 原 植被 NPP 
与 降水 .气温 的 相关 性 具有 一 定 的 空间 异 质 性 , 且 
二 者 的 相关 有 一 定 的 互补 性 , 即 与 降水 呈现 正 相 关 
的 地 区 和 和 气温 呈 负 相关 关系 ,与 降水 呈现 负 相 关 的 
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4 黄土 高 原 2000 一 2018 年 植被 NPP 与 降水 (a) 气温 (b) 的 相关 性 
Fig. 4 Correlation of vegetation NPP with precipitation (a) and air temperature (b) in the Loess Plateau from 2000 to 2018 


地 区 和 气温 呈正 相关 关系 。 
3.3 人 类 活动 对 植被 净 初 级 生产 力 的 影响 

植被 NPP 的 变化 是 气候 变化 和 人 类 活动 共同 
作用 的 结果 ,本 研究 基于 降水 .气温 两 个 气候 因子 ， 
首先 计算 植被 NPP 预 测 值 ,然后 用 植被 NPP 真 值 与 
预测 值 的 差 值 NPP 代 表 人 类 活动 对 植被 NPP 的 影 
响 。 由 图 $a 可 知 ,黄土 高 原 2000 一 2018 年 人 类 活动 
对 植被 NPP 的 影响 速率 每 年 为 3.72% ,呈现 由 负 转 
正 的 波动 上 升 趋势 ,表明 人 类 活动 对 植被 的 影响 由 
消极 变 为 积极 ,其 中 2000 年 和 2001 年 的 人 类 活动 
对 植被 NPP 主 要 表现 为 较 大 消极 影响 ,主要 是 因为 
退耕 还 林 还 草 工 程 实施 初期 ,退耕 地 上 没有 生长 出 
新 的 植被 , 即 人 类 活动 对 植被 表现 为 消极 影响 ;在 
2002 年 之 后 人 类 活动 对 植被 NPP 主要 表现 为 积极 
作用 ,是 因为 随 着 植被 的 逐渐 恢复 .环境 好 转 .生态 
得 到 改善 。 如 图 $b 可知, 人 类 活动 对 植被 NPP 影 响 


(a) 


4 1 1 1 1 4 4 1 J 
2000 2002 2004 2006 2008 2010 2012 2014 2016 2018 


年 份 


是 降低 的 趋势 面积 占 研究 区 的 15.75% , 分布 于 内 蒙 
古 北部 ,陕西 关中 地 区 和 河南 省 的 部 分 地 区 ;人 类 
活动 对 植被 NPP 影 响 表 现 为 增加 趋势 的 面积 占 比 
为 84.25% ,分 布 于 研究 区 的 大 部 分 地 区 ,其 中 陕西 
北部 的 延安 市 .榆林 市 地 区 人 类 活动 的 影响 最 为 明 
显 。 总 而 言 之 ,人 类 活动 在 黄土 高 原 植被 恢复 、 生 
态 环境 改善 的 过 程 中 具有 重要 影响 。 
3.4 气候 变化 和 人 类 活动 对 植被 净 初 级 生产 力 的 
相对 贡献 率 

黄土 高 原 的 植被 NPP 变 化 受到 气候 变化 和 人 
类 活动 的 共同 驱动 。 本 研究 基于 植被 NPP 真 值 、 预 
测 值 和 残 差 趋势 进行 相对 贡献 率 的 计算 ,得 到 
2000 一 2018 年 气候 变化 和 人 类 活动 对 黄土 高 原 植 
被 NPP 的 相对 贡献 率 分 别 为 48.789% .51.2296。 

如 图 6a 所 示 ,将 气候 变化 对 黄土 高 原 植 被 NPP 
的 贡献 率 分 为 5 个 等 级 :09%6~20% ,20% ~40% .40% ~ 


(b) 


Sz 


图 5 黄土 高 原 2000 一 2018 年 人 类 活动 对 植被 NPP 的 影响 
Fig. 5 Impact of human activities on NPP in Loess Plateau from 2000 to 2018 
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(a) 气候 变化 贡献 率 /% N  (b) 人 类 活动 贡献 率 /% N 
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6 2000—2018 年 气候 变化 (0 和 人 类 活动 (pb) 对 植被 NPP 贡献 率 的 空间 分 布 
Fig.6 Spatial distribution of climate change (a) and human activities (b) contributions to NPP from 2000 to 2018 


60% .60%~80% 和 80%~100% ,面积 占 比 由 少 到 多 为 
40% ~60% >20% ~40% >60% ~80% >0% ~20% >80% ~ 
100% ,气候 变化 贡献 率 在 60%~100% 的 占 比 为 
32.76% ,主要 分 布 在 河套 平原 一 带 、 山 西北 部 ( 朔 州 
市 大同 市) 和 南部 (长 治 市 .临汾 市 .晋城 市 .运城 
市 ) 地 区 ;其 次 气候 变化 贡献 率 在 20%~60% 的 面积 
占 比 最 大 ,主要 分 布 在 山西 中 部 (忻州 市 .太原 市 、 
晋 中 市 .阳泉 市 .吕梁 市 )、 陕 西北 部 (榆林 市 .延安 
市 )、 甘 肃 东部 (定西 市 临夏 市 )、 宁 夏 回 族 自治 区 
(固原 市 ) 和 青海 省 ( 海 东 市 ) ;气候 变化 贡献 率 为 
0%~20% 范 围 的 面积 占 比 为 17.22% ,主要 分 布 在 南 
部 庆 阳 市 .平凉 市 ` 宝 鸡 市 .吕梁 市 .临汾 市 的 部 分 
地 区 。 

如 图 6b 所 示 , 将 黄土 高 原 人 类 活动 对 植被 NPP 
的 贡献 率 分 为 5 个 等 级 :0%~20% .20% ~40% . 40% ~ 
60% .60%~80% 和 80%~100% ,面积 占 比 由 少 到 多 以 
HN 40% ~60% >60% ~80%>20% ~40%>80%~ 100%> 
0% ~20% ,其 中 人 类 活动 贡献 率 在 60%~100% 的 面 
积 占 比 为 38.85%, 主 要 分 布 在 山西 中 部 (如 梁 市 、 临 
汾 市 ) .陕西 (延安 市 .宝鸡 市 )、 庆 阳 市 .平凉 市 ;人 
类 活动 贡献 率 在 20%~60% 的 面积 占 比 为 45.89%， 
主要 分 布 在 山西 北部 和 上 毗邻 的 内 蒙古 自治 区 呼 和 
浩特 市 和 乌 兰 察 布 市 .甘肃 省 的 东南 部 、 陕 西北 部 ; 
人 类 活动 贡献 率 在 0%~20% 的 面积 占 比 为 15.26%， 
主要 零星 分 布 在 黄土 高 原 的 边缘 地 区 。 

在 省 域 尺度 上 ,气候 变化 对 植被 NPP 的 平均 贡 
献 率 为 47.98% ,贡献 率 由 小 到 大 为 甘肃 省 < 陕西 省 
< 河南 省 < 青海 省 < 山西 省 < 宁夏 回族 自治 区 < 内 


蒙古 自治 区 (图 7) ,其 中 只 有 甘肃 省 和 陕西 省 的 贡 
献 率 小 于 40% , 内 蒙古 自治 区 的 气候 变化 贡献 率 最 
高 为 64.4% , 其 次 为 宁夏 回族 自治 区 的 贡献 率 
53.83% 。 人 类 活动 对 植被 NPP 的 平均 贡献 率 为 
52.02% , 贡献 率 由 小 到 大 为 内 蒙古 自治 区 < 宁夏 回 
族 自治 区 < 山西 省 < 青海 省 < 河南 省 < 陕西 省 < 
甘肃 省 ,其 中 ,内 蒙古 自治 区 的 人 类 活动 的 贡献 率 
最 低 为 35.6% , 甘肃 省 的 人 类 活动 贡献 率 最 高 为 
60.79% ,其 次 为 陕西 省 的 人 类 活动 贡献 率 60.35%。 
气候 变化 在 人 口 稀 少 .发 展 相 对 落后 省 份 的 植被 
NPP 变 化 中 有 重要 影响 ;人 类 活动 在 人 口 众多 AGF 
发 展 较 好 省 份 植被 NPP 变 化 中 占 举足轻重 的 地 位 。 


E 气候 变化 mm 人 类 活动 


相对 贡献 率 /% 


陕西 ”山西 ”青海 ”宁夏 内 蒙古 河南 。 甘肃 
省 份 
图 7 黄土 高 原 2000 一 2018 年 不 同 省 份 气候 变化 和 人 类 活 
动 的 相对 贡献 率 
Fig.7 Relative contribution of climate change and human 


activities in different provinces 
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ar 水 土 流失 .保护 生态 环境 ,各 地 开始 实施 小 范围 的 
4 讨论 禁止 砍伐 放牧 、 大 规模 的 退耕 还 林 还 草 等 工程 3， 


4.1 黄土 高 原 植被 净 初 级 生产 力 时 空 分 布 格局 
成 因 

在 时 间 上 ,2000 一 2018 年 黄土 高 原 植被 NPP 的 
年 际 变化 较 大 (图 1) ,但 整体 上 呈现 波动 上 升 的 趋 
势 。 因 为 在 1999 年 后 实施 了 退耕 还 林 还 草 的 生态 
恢复 工程 ,而 黄土 高 原作 为 其 重点 示范 区 , 随 着 生 
态 工程 的 实施 ,植被 逐步 恢复 ,生态 环境 得 到 较 好 
的 改善 ,植被 NPP 呈现 上 升 的 趋势 ,这 与 Teng SE 
人 的 研究 结果 一 致 。 在 空间 上 ,2000 一 2018 年 黄土 
高 原 年 均 植 被 NPP 呈现 西 北 低 东 南 高 的 分 布 格局 
(图 2) ,黄土 高 原 86.86% 的 植被 NPP 旦 增加 趋势 , 主 
要 位 于 研究 区 中 部 的 陕西 北部 和 山西 吕梁 山区 , 即 
退耕 还 林 还 草 的 核心 区 域 ;13.14% 的 地 区 植被 NPP 
旺 下 降 趋势 ,主要 分 布 于 西北 部 的 干旱 沙漠 区 等 自 
然 环 境 较 差 的 地 区 ,这 与 Feng 等 中 和 Li 等 "的 研究 
结果 一 致 。 
42 黄土 高 原 气候 变化 及 人 类 活动 对 植被 净 初 级 
生产 力 的 作用 机 制 

气候 变化 会 使 植被 的 种 类 分布. 生长 状况 发 
生 改 变 以 适应 新 的 环境 站, 尤其 在 干旱 、 半 干旱 区 ， 
植被 对 气候 因子 变化 的 响应 更 为 敏感 。 如 图 3a 所 
AR ,研究 区 植被 NPP 与 降水 呈正 相关 关系 的 面积 
93.79% , 这 是 因为 黄土 高 原 是 典型 的 干旱 半 干 旱 
区 ,水 资源 是 限制 植被 生长 的 主要 气候 因素 ,降水 
的 增加 有 利于 植被 生长 汪 ; 呈 负 相关 关系 占 6.21%， 
主要 位 于 东南 部 的 汾 河流 域 .黄河 下 游 三 角 洲 地 
区 ,这 些 地 区 降雨 充沛 ,植被 变化 主要 受 人 类 活动 
的 作用 3。 如 图 3b 所 示 , 黄 土 高原 植 被 NPP 与 平均 
气温 的 相关 系数 为 -0.05 ,因为 温度 的 上 升 使 茸 散 量 
XE mn? 、 土 壤 水 分 不 足以 支持 植被 的 生长 ,甚至 导 
致 土壤 干 层 加 厚 ””。 植 被 变化 除 受 气候 变化 的 影 
响 外 ,还 与 人 类 活动 密 不 可 分 ”。 由 图 4 所 示 , 人 类 
活动 影响 呈 减 少 趋势 的 面积 占 比 为 1$% , 主要 位 于 
陕西 关中 平原 .内 蒙古 自治 区 北部 等 经 济 发 展 较 好 
的 、 人 口 众多 地 区 ,这 些 地 区 为 了 满足 城市 发 展 的 
需要 ,大 规模 的 耕地 等 被 转换 为 建设 用 地 ,导致 生 
态 环境 遭 到 破坏 ,人 类 活动 对 植被 NPP 表 现 为 消极 
作用 中; 人 类 活动 影响 呈 增 加 趋势 的 面积 
84.25% , 主要 位 于 陕 北 和 山西 吕梁 山区 ,为 了 缓解 


因此 , 陕 北 和 山西 吕梁 山区 人 类 活动 对 植被 NPP 表 
现 为 积极 影响 。 

本 研究 将 温度 .降水 作为 关键 气候 因子 ,基于 
残 差 分 析 剥 离 气 候 变化 和 人 类 活动 ,结果 显示 ， 
2000 一 2018 年 黄土 高 原 植被 NPP 变化 由 气候 变化 
(48.78% ) 和 人 类 活动 (51.22% ) 共 同 驱动 ,这 与 前 人 
的 研究 结果 一 致 上 和 。 此 外 , 残 差 分 析 虽 广泛 应 用 
于 分 离 气候 变化 和 人 类 活动 影响 因素 的 研究 ,但 其 
自身 也 存在 待 完 善 的 部 分 ,比如 ,本 人 研究 仅 选择 气 
温 、 降 水 代替 气候 因子 建立 与 植被 净 初 级 生产 力 的 
回归 方程 , 现 有 部 分 研究 ”I 也 考虑 了 太阳 辐射 等 
其 他 气候 因子 或 者 更 具有 生态 价值 的 气候 指数 等 ， 
未 来 研究 如 何 选择 气候 因子 增加 方程 的 合理 性 仍 
具有 不 确定 性 ;在 人 为 因素 方面 ,只 能 一 概 而 论 为 
人 类 活动 ,并 不 能 细 化 其 具体 方式 如 城市 扩张 、 生 
态 建 设 等 。 因 此 ,未 来 需要 通过 增加 气候 要 素 选 择 
的 合理 性 、 细 化 人 类 活动 的 具体 方式 ,达到 减少 植 
被 NPP 影 响 因素 研究 的 不 确定 性 。 


5 结论 


本 研究 基于 CASA 模型 估算 了 黄土 高 原 植 被 
NPP 数 据 集 ,通过 趋势 分 析 ` 相 关 性 和 残 差 分 析 等 方 
法 ,定量 探究 了 2000 一 2018 年 黄土 高 原 气 候 变 化 和 
人 类 活动 对 植被 NPP 的 影响 , 主要 结论 如 下 : 

(1) 黄土 高 原 2000 一 2018 年 植被 NPP 为 东南 
高 西北 低 的 分 布 格局 ,整体 上 呈 波 动 上 升 的 态势 ， 
研究 区 中 部 的 陕西 北部 和 山西 吕梁 山 呈 增加 趋势 ， 
西北 部 的 干旱 沙漠 区 等 自然 环境 较 差 的 地 区 植被 
NPP =Œ FBR, 

(2) 2000—2018 年 黄土 高 原 气 候 对 植被 NPP 变 
化 的 贡献 率 为 48.78% ,其 中 植被 NPP 与 降水 为 正 相 
关 关 系 与 气温 呈 负 相关 关系 ,气候 变化 起 主导 作 
用 的 区 域 主要 分 布 在 河套 平原 一 带 、 山 西北 部 和 南 
部 地 区 ;人 类 活动 对 植被 NPP 变化 的 贡献 率 为 
51.22% ,贡献 率 大 于 60% 的 区 域 主要 分 布 在 山西 中 
部 、 陕 西北 部 延安 市 和 中 部 宝鸡 市 、 庆 阳 市 . 平 
凉 市 。 

(3) 黄土 高 原作 为 干旱 半 干 旱 区 植被 对 气候 变 
化 的 响应 更 为 敏感 , 随 着 退耕 还 林 还 草 构成 的 大 规 
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模 实施 , 人 类 活动 在 生态 改善 方面 发 挥 着 重要 作 
用 。 本 文 仅 将 温度 、 降 水 作为 气候 变化 的 关键 因 
子 , 忽 视 了 其 他 可 能 对 植被 变化 产生 影响 的 因素 ， 
并 未 细 化 人 类 活动 的 具体 方法 。 因 此 ,黄土 高 原 气 
修 变 化 和 人 类 活动 对 植被 变化 的 影响 机 制 仍 需 次 
入 研究 ,以 期 为 区 域 生 态 恢 复 和 管理 提供 科学 
依据 。 
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Contribution of climatic change and human activities to changes in net primary 


productivity in the Loess Plateau 
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Engineering, Chang’ an University, Xi an 710064, Shaanxi, China; 3. Shaanxi Key Laboratory of Land Reclamation 
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Abstract: In this study, we used the CASA model to produce a net primary productivity (NPP) dataset for the 
Loess Plateau from 2000 to 2018. Based on correlation and residual analyses, we quantified the relative 
contribution of climate change and human activities to NPP, with a goal of providing a scientific basis for the 
formulation of local ecological protection strategies and improving watershed ecological security. The 
distribution pattern of NPP in the Loess Plateau from 2000 to 2018 was high in the southeast and low in the 
northwest. The portion of NPP that followed an increasing trend was mainly distributed in the core area of the 
Grain for Green Program, accounting for 86.86% of the total NPP, whereas the portion of the NPP that followed a 
decreasing trend was distributed in the northwestern arid zone, accounting for 13.14% of the total NPP. The 
contribution rates of climate change and human activities to NPP were 48.78% and 51.22%, respectively. Climate 
change plays an important role in vegetation change in sparsely populated and underdeveloped areas, whereas 
human activities play an important role in vegetation change in areas with large populations and better economic 
development. The Loess Plateau, an arid and semi-arid region, has a more sensitive vegetation response to climate 
change. With the increasing frequency of human activities, climate and human activities jointly determine the 
vegetation changes. This study contributes to the understanding of the effects of climate change and human 
activities on the dynamic changes of vegetation. It provides a scientific basis for vegetation restoration and high- 
quality development of the Loess Plateau. 
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